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We consider the phase-ordering kinetics of one-dimensional scalar systems. For attractive long- 
range (r~' 1+<T ') interactions with a > 0, "Energy- Scaling" arguments predict a growth-law of the 
average domain size L ~ f or a rj a > 0. Numerical results for a = 0.5, 1.0, and 1.5 

demonstrate both scaling and the predicted growth laws. For purely short-range interactions, an 
approach of Nagai and Kawasaki is asymptotically exact. For this case, the equal-time correlations 
scale, but the time-derivative correlations break scaling. The short-range solution also applies to 
systems with long-range interactions when a — » 00, and in that limit the amplitude of the growth 



law is exactly calculated. 
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I. INTRODUCTION 

> 

^-j- ■ The quench of a system from a disordered phase to an ordered phase is a non-equilibrium process in which energy 
is dissipated and any topological defects are eliminated. Typically, the system develops a scaling structure, with a 
single time-dependent length scale, as the growing broken-symmetry phases compete to select the ordered phase [[l]. 
In this paper we discuss non-conserved one-dimensional (ID) scalar systems, where the competing phases are domains 
of uniform magnetization and the topological defects are domain walls. We limit our discussion to the late stages of 
phase ordering, when the domains are much larger than the interface width and the scaling structure, if any, has been 
established. The phase-ordering dynamics is driven by interactions between domains, rather than by the curvature of 
domain walls which drives the dynamics of short-range systems in higher dimensions [g| . For short-range ID systems 
this interaction is through the exponential tail of the domain- wall profile [^). For long-range systems, the domains 
have a direct r~( 1+CT ' interaction. 

Following work in dimensions greater than one , there has been recent interest in ID systems with attractive 
long-range interactions. This paper follows up a broader treatment that determined growth laws for for a wide variety 
of scaling systems, including ID systems with < a < 1 |Q. Lee and Cardy || found analytic and numerical 
^ , indications of scaling violations and anomalous growth laws for < a < 1, in addition to results consistent with 
scaling for a > 1. This motivated us to extend their numerical studies for systems with a — 0.5, 1.0, and 1.5, and to 
extend our theoretical treatment to all a > (section |J). We find that larger systems than those studied by Lee and 



Cardy exhibit both scaling and our predicted growth laws (section II 13). 



Short-range ID scalar systems have been studied both numericall y |S| j and experimentally |10[ . Nagai and Kawasaki 
have presented a solution assuming uncorrelated domain sizes [ p"T| , |12[ . We rederive their solution with a simplified 
approach, and show that it is asymptotically exact for systems with purely short-range interactions, and also for 



systems with long-range interactions in the limit a — > 00 (section HI) 



A generic energy functional for a ID system with a scalar order parameter, <j>(x), and short-range interactions, is 

H[(f>}= I ' dx[{d(j>/dxf + V(cf>)] , (1.1) 
where V(<f>) is a double well potential such as 

V(4>) = (0 2 - l) 2 . (1.2) 

We restrict our attention to a potential with symmetric quadratic minima Jl3[ . We can also add long-range attractive 
interactions, 
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(1.3) 

where is the system size and we take a > for a well defined thermodynamic limit. The long-range interactions 
dominate the ID dynamics when they exist, so we have ignored the short-range component of (1.3) for a > 2. 



After a temperature quench into the ordered phase, the equation of motion for the ordering kinetics for non-conserved 



systems with purely dissipative dynamics is 



opt 




d t 4> = -5H/S4>. 



(1.4) 



For short-range systems, and long-range systems with a > 1 |lq| , there is an ordered phase only at T = 0. We restrict 
ourselves to T — 0, with no thermal noise, for all cases, since temperature is an "irrelevant variable" within an ordered 
phase |16|. 



II. LONG-RANGE ATTRACTIVE INTERACTIONS 



We first treat systems with long-range attractive interactions of the form (1.3). As was shown by Lee and Cardy 
for Ising spins B, the dynamics of the system is captured by the motion of the domain walls. Consider two sharp 
domain walls of the same sign at Ti and rj , as shown in figure [l| Their interaction energy, Eij , is the part of i?LR 
that depends on their separation: 
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(2.1) 



where we have scaled out a factor of 8/cr. Domains of opposite sign introduce an overall —1 factor to the energy. This 
interaction energy leads to a force, 



F tj = -dEij/drj, 



(2.2) 



where I 



Ti\ and the force is repulsive between walls of the same sign and attractive between walls of the 



opposite sign. The domain wall velocity is proportional to the summed pairwise forces acting on it [Jlj 

n = EX 



(2.3) 



When adjacent domain walls meet, they annihilate — this drives the non-equilibrium coarsening process. We use the 
domain-wall dynamics described by equations ( [2.2] ) and (2.3) unless spec ifica lly noted. 

This domain-wall dynamics also hold s with the soft potential V(tfi) in ( |l.2| ). Fixing the domain wall positions {r^}, 
we use the quasi-static solution of (1.4) expressed in real space: 
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(2.4) 



where x 6 (r^r^+i). We use the saturated values <p rj ±1 within the integral, so that the integrals are restricted to 
phases with 4>(x') = —<f>(x). We also assume that x is far from all domain walls, so that we neglect dl<fi. This leads 
to <f>(x) = ±(1 — p(x)) with 
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(2.5) 
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at large distances from other domain walls (compare M). The force corrections due to the asymptotic profile can 
be easily estimated by the energy per unit length of a dding to the ends of the power-law profile between two defects 



as they move apart. From the potential term in (1.1), V(<fi) = (<fi 2 — l) 2 ~ l~ 2a , where the defects are a distance I 
apart. From the gradient squared-term, the effective force is l~ 2 ~ 2rT , which at late times is always much less than from 
the local potential. Since we are assuming a quasi-static configuration, the change in energy due to the long-range 
interactions is of the order of that due to the potential term, so is also l~ 2a ' . These corrections are negligible in 
comparison to the l~ a interactions of "sharp" domain walls (|]^), which we use without loss of generality. 



A. "Energy-scaling" theory for growth laws 

We predict the growth-law of the characteristic length, L(t), for scaling systems by considering the energetics of 
the system 0. The sca ling behavior, or dependence on the length-scale, of the energy density is described by the 



long-range interactions (1.3) 



dkk a (faf-k) , (2.6) 

where the angle brackets indicate an average over initial conditions. We independently calculate t he r ate of energy- 
density dissipation by integrating the dissipation in each Fourier mode and using the dynamics ( JxT3 ) expressed in 
momentum space: 



de/dt = I dk ((SH/S^ k )d t (t> k ) , 

dk {d t 4> k d t 4>- k ) . (2.7) 



We use the equal-time correlation function, S(k, t) = (<fik<fi-k}, to determine the energy density, and the time-derivative 
correlation function, T{k,t) = (dt4>kdt4>-k) , to determine the rate of energy-density dissipation. 

If we assume scaling — that the correlations are described by one time-dependent length-scale L(t) but not by any 
microscopic scales — then dimensional analysis determines the structure function: 

S(k,t) = L(t)g(kL(t)), (2.8) 

where the scaling function g(x) has no time-dependence. This is also known as the dynamic scaling hypothesis 
p8| . We apply the same scaling hypothesis to the two-time correlation function: S(k,t,t') = {<j)k{t)(j)-k{t')) = 
&T~\g(fcL(i), kL(t'), t/t'). This leads to a scaling form for the time-derivative correlation function: 

(d t <p k d t ^ k ) = Wt=t'S(fc,M0, 

= (L/L) 2 Lg xy (kL) + r 2 Lg z {kL), (2.9) 
~ (L/L) 2 Lg tt (kL), (2.10) 

where for power-law grow th, with o r wi thout logarithmic factors, the g z contribution does not dominate at late times. 



If the energy integrals (2.6) and (2.7) are convergent then we use the scaling forms and change variables to get the 
L and L dependence. Certainly, the integrals are convergent for kL -C 1 if the thermodynamic limit exists. However, 
the integrals may diverge for kL ^> 1. For S(k,t), that limit is given by Porod's law |19| , in which the structure at 
kL ^> 1 is proportional to the wall density p wa \\ ~ L , 

S(k, t) ~ L" 1 k~ 2 , fcL>l. (2.11) 



With this and the scaling form (|2.8|) we determine the energy density (2.6): 




a > 1, 

<7=1, (2.12) 
a < 1. 

For a > 1 the integral diverges for kL 3> 1 and the energy density is proportional to the domain-wall density, 
Pwaii ~ L~ , for cr = l the logarithmic divergence of the integral indicates that structure on all scales contributes to 
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the energy-density; while for a < 1 the integral converges and the energy is dominated by the long-range interaction 
between domains. 

The interaction energy-density, e; nt , of domains at the characteristic scale is the convergent kL rj 1 part of the 
energy integral. This energy-density of interaction is what drives the motion of well separated domain walls: 

e lnt ~ L- a . (2.13) 

To determine the growth law, we only need the time-derivative of e- m t, which comes from domain walls moving at the 
characteristic velocity L. Near a particular domain wall the order parameter comoves with the domain wall: 

d t <j> = -vd<t>/dx. (2.14) 

The contribution of these slowly moving domain walls in momentum space is 

{d t 4> k d t ^ k ) int ~ L 2 k 2 <^_ fe ) , 

~L 2 /L, fcL>l, (2.15) 



where we use Porod's law to obtain the second line. While this expression satisfies the scaling form (2.10), it is not the 



full time-derivative correlation fu ncti on, merely the part that contributes to the dissipation of the interaction energy 



£int- We use this contribution in ( |2.7| ) to find 

de int /dt ~ r 1 L 2 /L. (2.16) 

The integral diverges for kL ^ 1, but the divergence reflects the internal structure of the domain walls, rather than 
small separations with respect to L(t). We compare the rate of energy-density dissipation to the time derivative of 



(2.13) to obtain L ~ L <T , which determines the growth law: 

L(t) (2.17) 

for all a > 0. 

We check the consistency of this approach by calculating the full de/dt. This is easily done in real space, where 
from (2.7) we have 

de/dt = -L^ J dx {dt(f){x)dt(j){x)) , 

~ L^ 1 (y 2 ) , (2.18) 



where we have used (2.14) to obtain the second line. The result is the density of domain walls, /0 W aii ~ L 1 , times the 
average square-velocity. Either (u 2 ) = L 2 and the energy-density is the same as ei n t, or the energy is dominated by 
the domain wall energy and the energy-dissipation is dominated by rapidly annihilating domain walls. 

For domain sizes / <C L, the force due to nearest-neighbors will set the scale of the net force on a domain wall, so 
from ( [2~2| ): 

K0| ra (3l- a , Z<L, (2.19) 

where the constant prefactor (3, where < (3 < 1, proves sufficient to parameterize any shielding due to other domain 
walls. If n t (l) is the number distribution of domain sizes, so that n t (l)Sl is the number of domains of sizes in the 
interval [1,1 + SI], then the number flux of collapsing domains is given by = n t (l) [— 2\v(l)\] for I <§; L. The factor of 
two follows since the collapse rate of a small domain is twice the speed of each of its walls. When a domain collapses 
its domain walls vanish, merging the two adjacent domains. Hence, every domain that collapses decreases the total 
number of domains by 2. Matching 2j(£) to the rate of change of domain number gives 2j(£,t) = dtN = —LL^jL 2 , 
where N = L^jL = J, dlnt(l). Since domain walls annihilate at a constant rate over timescales much less than 

L/L, the domain wall flux is constant for scales I -€i L: « j(£). This leads to 

N ■ 

rn^K—Ll*, /«L, (2.20) 

where /3 = 1 i f no s hieldi ng of small collapsing domains occurs. 

With ( 2.19| ) and (2.20) we determine the contribution of domains of scale I <C L to the average square velocity: 
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j L dlv 2 (l)n t (l) 
(v ) n = — — f 

\ / small rL / 

L dln t (l) 




a > I, 

a=l, (2.21) 
a< 1. 

For cr < 1 small domains do not dominate the average and, since L ~ we have (u 2 ) ~ L 2 . Indeed, for this 

case e ~ £i n t, and the energe tics i s contr olled by domains at the c hara cteristic scale. For a > 1, the energy-density 
dissipation, calculated with ( [2.1§| ) and ( 2.21 ), is consistent with ( 2.12j ). However, the energy-density dissipation is 



dominated by small annihilating domains and does not determine the growth law of the average domain size, for 
which we have used e;„t. 

We have numerically studied the average square velocity of domain walls, as well as the average velocity. In Figure 
H we show (v 2 } / L 2 and (\v\) /L vs. time. The scaled velocity is approximately constant away from early time 
transients and late time small number effects. Indeed, a similar analysis to that above shows (\v\) ~ L for all a > 0. 
The average square-velocity scales well for a = 0.5. However for a = 1.0 and 1.5, (v 2 ) is much larger than L 2 - 
corresponding to the dominance of small rapidly annihilating domains for a > 1. While sparse statistics for extremely 
small domains leads to sporadic undercounting of (i> 2 ) for a > 1, the scaled square-velocity for a — 1.5 agrees with 
(y 2 ) /L 2 ~ (LL)^ 1 ~ i 1 / 3 , which is shown by the straight line. 



B. Simulations 



To check these growth laws and the scaling of the correlation functions, we perform computer simulations for 
a = 0.5, 1, and 1.5. These are the cases considered by Lee and Cardy and our treatment is motivated by theirs. 
We simulate the evolution of systems with sharp domain walls of alternating sign which annihilate when they meet. 
These domain walls are treated as particles. They interact with long-range forces, f(l) = H~ a , where I is the domain 
wall separation. The force is repulsive or attractive for like or unlike domain walls, respectively. 

We first measure the number of domain walls N(t) as a function of time, starting from a fixed numbe r pl aced 
randomly |20| , and continuing until none remain. With a domain wall velocity proportional to the net force ( |2.3| ), we 
use a simple Euler update for each domain wall location, Ar^ = r.;Ai. Only the early-time behavior depends on the 
size and nature of the timestep At, if At is small enough. We use a timestep proportional to the time (At — i/20), 
which for power-law growth results in a fixed minimum scaled domain size that annihilates in one timestep pjj], 

U„/Loc (AiA) 1/(1+ff) . (2.22) 

The choice of boundary condition only affects the late-time results. We use free boundary conditions |2^] for which 
small number effects are not seen until the average number of domain walls N < 1. We average over at least a 
thousand independent runs for most systems, and over at least 50 runs for the largest system size with 12800 initial 
domain walls. 

We expect the average domain size L(t), which is the inverse of the domain wall density pit), to have our predicted 
growth law L(t) = £ 1 /( 1 + <T ) /<5>*(er). The ratio of the measured domain wall density to the predicted, p(t)t 1 ^ 1+a ' ) up 
to the amplitude factor, should be a scaling function of dimensionless combinations of lengths: 

p(t)tVd+-) = $ [ p (t)/p , l/ipm^)} . (2.23) 

In Figure [| we plot, for various a , the scaling function $ as a function of the number of domain walls, N, on the left 
and of p(t)/p on the right. We see an early-time transient (i.e. po/p(t) of order 1) which is independent of the system 
size, and a small number regime (i.e. N = p(£)£oo small) at late times which is independent of the initial density [po[ . 
At intermediate times the scaling function is constant, $ = Q*(<r) = $[0,0], independent of both the initial density 
and the system size. This indicates a scaling regime, growing with system size, that confirms our predicted growth 
law L(t) = < 1// ( 1+CT '/<f>*((T). In table | we plot the amplitude Q*(a), and also the fixed point amplitude used by Lee 
and Cardy g* R (a) = $*(cr) [2(1 + cr)] l7(1+CT) §. We include the exact amplitude for cr — > oo that we derive in section 
III (see equation |3.15| ) . 
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In Figure [|, we plot $ vs. N with two different boundary conditions and a = 0.5. The scaling regime is of different 
size in the two systems, but the scaling value $*(cr) is unchanged. The difference is only in the small number regime 
(N < 20). We note that systems with periodic boundary conditions and less than 800 initial domain walls, as studied 
by Lee and Cardy for a — 0.5 [0, do not probe the scaling regime. We find the expected scaling growth law, L ~ t 2 ^ 3 
for a = 0.5, in contrast to their numerical results. 

The ratio of the measured energy-density, e(t), to the predicted, e (2.12), should be a new scaling function of 
dimensionless combinations of variables: 

e(t)/e = <& B [p(t)/p , l/(p(t)£oo)] • (2-24) 

In the scaling regime we expect $e = $b[0, 0] to be constant. We plot e(t)/e vs. p(t)/p in Fig || for a = 0.5 
and 1.0, where we calculate the predicted e, up to a constant factor, from the measured length-scale L(t) = L^/N. 
We use periodic boundary conditions to avoid strong end-effects in the energy. The domain walls interact with an 
infinite-ranged force law so that all images of a domain wall around the periodic system are seen (see Appendix |A|). 
Note that only 21 runs are done for the largest a = 1.0 system. We see the expected early-time transient independent 
of system size, an intermediate scaling regime growing with system size, and a small-number late-time regime. The 
energy density for a — 1.5 is dominated by the core energy and is simply proportional to the number of domain walls. 

While the growth laws agree with the scaling predictions, to investigate scaling we must consider the scaling of the 
correlation functions directly. If scaling holds, the "Energy-Scaling" argument determines the growth laws, but the 
growth laws are not sufficient to determine scaling. To study correlations, we use periodic boundary conditions with 
infinite r ang ed forces. We also use initial conditions randomly chosen from the a — » oo fixed-point distribution (see 
equation 3.9), which reduces the initial transients for all a investigated. We generally use At = t/20, and average over 
at least 10000 independent runs, but we present some data at small scaled distances with a finer timestep, At = i/80, 
where we average over 1000 runs. 

We first plot the scaled number distribution of domain sizes, Fl(1/L) = Ln t (l)/N, against the scaled domain size, 
x = l/L, in Figure EL Th e exa ct result for a — > oo is also shown. The straight lines are the fit, for small x, to the 
power-law behavior from ( |2.20|) . The deviation from the power la w at very small x is a result of the geometrically 
increasing time-step in the simulations. Domains of size below i m ; n (2.22) tend to be randomly distributed. The filled 
points have four times as many updates per timescale, and are correspondingly closer to the power-law. We tabulate 
the fitted values of (3 in table O. The decreasing values of (3 for decreasing a indicate stronger shielding and greater 
correlations between domains. 

To obtain the correlation functions, we relate the gradient of the field to the positions of the sharp domain walls: 



JY 



d x cl>(r) = 2j2S(r-r i )(-iy, 



(2.25) 



where the sum is over all domain walls, of alternating sign, at positions {r{\. This equation gives the correct 
discontinuity of the field, 2(— l) 1 , at the domain walls. Similarly, the time derivative of the field is 



JV 



(2.26) 



where i?j is the signed velocity of the domain wall. We then easily calculate 
C"(r,t) = -(d r ct>(r + r )d ro <i>(ro)} rn , 



= -4L: 



ro 



ri)S(r 



and 



-4S(x) - 2N- 1 ^(-l) l+3 {S(x - \x i:j \) + S(x + \ Xij \)} 



(2.27) 



T(r,t) = (d t 0(r + r o )d t <f>(ro)) ro , 

= 4L- 1 J dr ^(-l) 4+ ^<5(r + r - n)6(r - rj ), 



(L/L)* 



4 (v 2 ) L- 2 5{x) + 2N- 1 Y^{-^Y +3 v l v 3 L- 2 {S(x - \x tJ \) + S(x + \ Xij \)} 



(2.28) 



G 



where x = r/L, Xij = (rj — ri)/L, — NL(t), and (v 2 ^ is the average square velocity of a domain wall. The 
real-space scaling forms are 

C"(r,t)=L- 2 f"(r/L(t)), (2.29) 



T(r,t) = (L/L) 2 g(r/L), 



(2.30) 



where the scaling functions f"(x) and g(x) are time-independent if scaling holds. Wc plot f"(x) in the top row and 
g(x) in the bottom row of figure for a = 0.5, 1.0, an d 1.5 . Shown are two times differing by a factor of 4, both 
chosen from the scaling region of the scaling function $ ( |2.23 ). There are no discernable scaling violations. The solid 
curves through the data are parameter less estimates assuming independently collapsing domains, which should be 
valid for small x > 0: 



f'(x)^2N- 1 J2 S(x-\xij\), 



<ij> 



±F L (x), 



(2.31) 



and 



g(x) » 27V- 1 [vij/i] 2 S(x- \ Xij \) 



v(xL)/i] F L {x), 
4/3 2 [^(a)f 1+(T) (l + a) 2 x- 2 -F L (x), 



x < 1, 



(2.32) 



where the sum is ov er ne ighboring domain walls, $*(cr) is measured from figure [?| and Fl(x) is taken from figure]^. 
Following equation ( [2.19] ), we have used v(l) — (3l~ a with (3 from table We have also used the scaling amplitude of 
the growth law, L(t) = t 1 ^ 1 ^ /$*(a). For f"( x), t he estimate fl2.3l| ), assuming independently collapsing domains, 
is excellent at small x. For g(x), the estimate ( 2.32 ) agrees at very small x. An intermediate regime, which grows 
with decreasing a, deviates from the estimate and indicates growing corre lation s between domain wall velocities in 
different domains. Note that due to the geometrically increasing timestep ([2.22 ) in our simulations, the domain size 
distribution deviates from the expected Fl(x) ~ x a for i C 1 (see figure H). For At/t — > 0, f"(x) and g(x) s hould 



follow the power-laws indicated by the offset straight lines in the same regime where the estimates ( 2.31 ) and ( 2.32 ) 
hold. This is evident for f"(x) with the data for At/t = 1/80, but a fixed time-step is needed to check the indicated 
power-law regime in g(x). 

We have not shown the 4 (w 2 ) L~ 2 5(x) contribution to g(x). From equation (2.21), we see that this breaks scaling 
for (j > 1 due to small rapidly annihilating domains (see also figure ^). Hence, T(k,t) has a constant contribution 
that breaks scaling for all k for a > 1. However, this does not change our "Energy-scaling" argument, which depends, 
through ei n t, only on domains of the characteristic size. 



III. SHORT-RANGE INTERACTIONS: CONTINUOUS SPIN MODEL 



For systems with only short-range interactions, the domain wall profiles lead to interactions between adjacent 
domain walls which drive the phase-ordering. From the dynamical equation ( |1.4| ), the equilibrium profile of an 
isolated domain wall is determined by d 2 (f>/dx 2 = 40(0 2 — 1) — which gives an exponential tail to the domain wall 
profile far from the domain wall. The domain walls, treated as particles, are attracted by a corresponding exponential 
interaction V(r) — Voe~\ r \^, where £ is the effective width of the domain wall p), and a resulting pairwise force 
F(r) — —dV/dr = Vo/£e-' r '/'». The domain wall velocities are then simply proportional to their net force. This 
case has been studied numerically [9J, and has been solved by Nagai and Kawasaki Jll| with the assumption of no 
correlations between domain sizes in the fixed point distribution. We rederive their results using a simpler approach 
p3fl . We also show that the assumption of uncorrelated domain sizes is true at late times, so that the solution of 
Nagai and Kawasaki is exact. Their solution also applies to systems with long-range interactions, in the limit a — > oo. 

When a domain annihilates, the two adjacent domains combine into one large domain which is the size of the 
original three domains put together. For systems with short-range interactions, or for long-range interactions with 
a — > oo, the shortest domain collapses instantaneously with respect to longer domains (see Figure^). For short-range 
models, size differences between nearby domains are of the order of the growing length-scale L(t) and shorter domains 
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collapse exponentially faster than any larger ones, with an exponent equal to the size difference divided by the fixed 
core-scale £. For long-range mode ls, the smallest domain annihilates instantaneously relative to longer domains for 
cr — > oo due to the pairwise force fl2.2|). 

Consider a random initial set of domain lengths. The collapse of the shortest domain corresponds to taking the 
smallest length from the set, adding it to the lengths of its neighboring domains, and replacing the three lengths by 
one length equal to their sum. As this procedure is iterated, the distribution of domain sizes, scaled by the size of 
the smallest domain, L(t), will approach a fixed point distribution. There will be no correlations formed between 
the lengths of different domains because the collapse of the shortest domain forms a single domain with its neighbors 
independently of the other domains in the system. This is similar, in that respect, to the "paste-all" model of Derrida 
et al. [p3| . Since all higher-point correlations are zero, the domain size distribution function describes the fixed point. 

If n t (l)Sl is the number of domains of sizes in the interval [/, I + 81] in the entire system, then 

/>oo 

N(t) = J n t (l)dl, (3.1) 



Loo = y_ ln t (l)dl, (3.2) 

f t (l)=n t (l)/N(t), (3.3) 

F(l/L,t)=Lf t (l), (3.4) 

where N(t) is the number of domains in the system, L M is the length of the system, f t (l) is the probability density 
of domain sizes, and F(x,t) is the s cale d probability density with respect to the minimum domain size L. Note that 



nt (I) is the same as used in section [I A . We will determine the time- independent fixed point distribution F(x). 

Each domain that annihilates combines with the two adjacent domains, each of length I with probability density 
ft(l)- Conversely, for each domain that annihilates, one large domain is formed at length I with a probability density 
coming from two domains combining, one of length I' and the other of length I — I' — L, both restricted to be larger 
than L. In a time interval 5t, the minimum length will shift by SL = L(t + St) — L(t), and nt(L)8L domains will 
annihilate. This leads to an evolution equation for the number distribution n t {l): 



n t +st(l) = n t {l) + n t (L)SL 



(3.5) 



-2/t(i) + y. dl'f t (l')f t (l-l'-L)Q{l-l'-2L) . 
The evolution equation for the scaled probability density F(x,i) follows: 

LdF(x,t)/dL = F(x,t) + xdF(x,t)/dx + F(l,t) / dx'F(x', t)F(x - x' - 1, t)Q(x - x' - 2). (3.6) 



Using the Laplace transform, (j)(p,t) = dxe px F(x,t), we obtain 

Ld<p/dL = -pd(f)/dp + F(l, t)e~ p [<p 2 - l] . (3.7) 
For the stationary solution we set d<p/dL = 0, use 0(0) = 1, and force consistency on d<f>/dp\ p= o to determine 



F(l) = 1/2 ||24| . Hence the stationary solution to ( p.7| ) is 



0(p)=tanh[Ei(p)/2], (3.8) 

where Ei(p) = dxer x jx. This agrees with Nagai and Kawasaki p| . Expanding the tanh and doing the inverse 
transform leads to the fixed point distribution 

[0b-1)/2] .ooSi+l , _^ 

F(x)= J2 IT— *(*-!>)' a^ 1 ' ( 3 - 9 ) 

i=0 Jl 3 = 1 Xj 

where the first sum is for all i such that x > 2i + 1, and 2 n A n are the Taylor expansion coefficients of tanh(x). The 
fixed point distribution is piecewise analytic, with discontinuities at x — 2i + 1 in the 2i-th derivative for all integer 
i > 0. The fixed point distribution is simple for small x: 



8 



(3.10) 



but gets progressively more complicated in every succeeding interval of width 2. 

At finite times, for systems with short-range interactions, there will be small correlations between the sizes of 
different domains because the exponential suppression of the motion of domain walls relative to those of the smallest 
domain is not complete. However, the ratio R of the average force on a domain wall to the force on the walls of the 
shortest domain vanishes as L(t) — > oo, and so any correlations will vanish in the scaling, or late time, limit. For 
short-range forces, 

R~ / dxF{x)e- Lx/ Z/e- L/ t 7 



1 



/i 2x 

£/(2L), (3.11) 



where F(x) is the fixed point distribution of the scaled domain size (3.9). Hence R vanishes at late times as L — > oo, 



and the fixed point distribution will hold in that limit. For long-range forces, the same calculation leads to 

/oo 
dxF(x)(Lx)- a /L- a , 



OO T_ 

dx 
i 2x* 



1/(2<t). (3.12) 



In the limit a — > oo this ratio vanishes and the fixed point distribution is valid at all times. 
The average domain size is given by 



/oo 
xF(x)dx, 



= -Ld(j)/dp\ p= o, 

= Z2e 7E , (3.13) 



where we have used the asymptotic form of (3.S), and — 0.577 is the Euler constant. This agrees with Nagai and 
Kawasaki |T]J and numerical results ||. 

The scaled distribution of domain sizes, Fl{x) — Ln t (xL)/N, is plotted against the scaled domain size, x — l/L(t), 
for a — 0.5, 1.0, and 1.5 in Figure |[ The exact result for a — > oo is also shown. As a increases, the distribution 
approaches the a — > oo fixed point distribution. 

The growth law of the average domain size, L(t) = L^/N^), is derived by matching the flux of annihilating domains 



j(£) to the flux of domains coming from the minimal domain size L. From the discussion leading to equation (2.20), 
the flux of annihilating domains is given by = N/2 = —LN/(2L). For systems with short-range interactions, the 
flux at L is given by j(L) = — (N/L)e~ L ^ , where n t (L) = (N/L)F(1) is the number density of domains at scale L, 
they are collapsing at a rate —2e^ L ^ (twice the speed of each end), and -F(l) = 1/2. Equating these fluxes give a 
logarithmic growth law 

L = £ln(t/to), (3.14) 

where to depends on the form of the potential V ((/))■ 

For systems with long-range interactions, with a — > oo, we can also obtain the exact growth law. The flux 
of annihilating domains is again given by = —LN/(2L). The flux at L is now given by j(L) = —N/LL a , 
where the domains of scale L are collapsing at a rate —2L~ a . Equating these fluxes gives dL/dt = 2L~ a , and so 
L(t) = 2e TB i 1 /( 1+cr ) in the a — > oo limit. This determines the scaling amplitude in that limit: 

$*(oo)= 5 ^(oo)=e- 7B /2~0.28, (3.15) 

where the fixed point amplitude used by Lee and Cardy is g* R {a) = $*(a) [2(1 + a)] 1 ^ 1+r7 " 1 M. This exact result is 
comparable to the perturbative calculation by Lee and Cardy who found <?|j(oo) « 0.33. 



9 



The structure function is calculated by Kawasaki et al. |L2] from the Fourier transform of (2.27) and the domain 
size distribution function: 

2 



S(k) = 



1 



tanh 



Ei(ikL)/2 



k 2 L 



1 + tanh 



Ei(ifcL)/2 



(3.16) 



which, because L/L is constant, satisfies the scaling form (2 



Since domains annihilate indepen dently and only adjacent domain walls have correlated velocities, the real-space 
time-derivative correlation function ( 2.28 ) is simplified 



T SR (r,i) = (L/L) 2 \4(v 2 )L- 2 S(x)+4N- 1 ]T (-ly^v^L^SQx] - \ Xij \) \ , 



<i]> 



(L/L) 2 4 (v 2 ) L- 2 S(x) + 4F L (\x\) v(\x\L)/L 



(3.17) 



where x = r/L, Fl(x) is the domain distribution function scaled with respect to the average domain size L, and 
(v 2 ') is the average square velocity of a domain wall. The equation becomes exact, though the scaling behavior is 
unchanged, if we use the RMS velocity of a wall of a domain of size I in the second line, wrms = |v(0 — (v) |> rather 
than v(l). For both short-range and a — > oo s ystems, the average square-velocity is dominated by the tiny flux of 
annihilating domains. Following equation (2.21), the average square- velocity will be 

(v 2 ) ~ LL- 1 dlv(l). (3.18) 

The integral diverges at small I for both short-range systems and long-range systems with a > 1, and the square- 
velocity is dominated by small annihilating domains, (u 2 ) ~ L/L. Hence, the scaling ( 2.10 ) of the time-derivative 
correlations (3.17) is broken at x = 0. The scaling for x > holds if v(xL)/L is time-independent. For long-range 
interactions, with v(l) ~ l~ a and L ~ L~" , this is true and the T(r, t) scales for r/L > 0. For short-range interactions, 

-l/xi 



V(l) 



Since L = £ln (t/to) (p-14|) , scaling only holds at x = L/L ~ 0.28, and is broken at all other x. 



IV. CONCLUSIONS 



From the energetics of scaling ID scalar systems, we predict growth laws L(t) ~ 



growth law is of the same form as the length-scale given by a collapsing isolated domain (2.2 



for all a > 0. This 
The intuitive picture 

is of collapsing domains leaving "voids" which set the growing length-scale. These growth laws agree with previous 
predictions for < a < 1 [^j]. We confirm our predicted growth laws, with simulations for a = 0.5, 1.0, and 1.5 
— one value in each of the three regimes of the equilibrium system. Different boundary conditions and system sizes 
only affect the late time, small number regime of our simulations, and so the thermodynamic limit is well behaved. 
The anomalous growth law seen by Lee and Cardy || for a = 0.5 was the result of large finite-size effects — a 
scaling regime appears for larger systems than they considered. The non-analyticities seen by Lee and Cardy || 
for < a < 1 must be controlled in an infinite system with a fixed initial density. We find that the equal-time 
correlations scale for a = 0.5, 1.0, and 1.5, as do the time-derivative correlations for a — 0.5. For a > 1 the real-space 
time-derivative correlations scales properly for r > £, however scaling is technically broken at core-scales due to small 
rapidly annihilating domains. For £ < r « L, the correlations can be well approximated by assuming that domains 
collapse independently. The approximation becomes better with increasing a. 

For purely short-range interactions, and for long-range interactions in the a — > oo limit, we show that the solution 
by Nagai and Kawasaki [ pd| of the domain-size distribution, which assumes uncorrelated domain sizes, is exact. Hence 
their expression for the structure factor ( 3.1 6[ ) , which satisfies scaling, is also exact. We determine the time-derivative 
correlation function, which breaks scaling at core scales, and more gen erally for systems with purely short-range 
interactions. For systems with long-range interactions described by (2.2) and (2.3), we determine the exact growth 
law amplitude in the a -> oo limit, L(t) = t 1/(1+ °' ) /fl , H (°°) with 9rI°°) = e ~ 7E / 2 - °- 28 - 
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APPENDIX A: RESUMMING DOMAIN- WALL INTERACTIONS IN A PERIODIC SYSTEM 



With long-range forces in a periodic sys tem , each domain wall interacts with an infinite number of images of each 
other domain wall. Using the interaction (2.2), the total force between two domain walls is given by 



/(0 - + J2 [( nL oo - iy a - (nioo + o _CT ] , 



n=l 
oo 



= -r- + ^[r( ( r)]- 1 

= -r a 



n=l 
t —a r°° 



dtt c 



fa ,-l^Mf/M e - I /2 ] 

sinh [a: / 2] 



(Al) 



where I is the closest distance between the domain walls around the loop and is the system size. The force is 
attractive along their nearest separation if the domain walls are of opposite sign and repulsive if of the same sign. 
The interaction energy can be similarly calculated using (2.1): 



E{l)-E(L 00 ) = l 1 -° + 



Ll 



IV - 1) Jo 



dxx 



cr-1 



cosh(/x/L c 



1 



sinh(a:/2) 



-a/2 



(A2) 



where we have adsorbed a factor of (a — 1) . An overall —1 factor is introduced between domain walls of opposite 
sign. For u / 1 the integrals are done numerically and stored in lookup tables for use in the simulations. For a = 1, 
the integrals can be done exactly: 



fa=l{l) = -7Ti 00 1 / tan ( 7r V L oo), 



(A3) 



E <r=1 (l)-E a=1 (L 00 )=lnl + hx 



sm(irl/L c 



wl/L c 



(A4) 
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FIG. 1. A schematic representation of the ID scalar system with domains of <f> = ±1 shown. The domain walls at n and 
Tj are "positive", and the sign of domain walls alternate along the system. Domain walls of the same sign repel each other, 
while those of opposite sign attract each other. The dynamics for long-range interactions are described by these sharp Ising-like 
walls, as is discussed in text. 



FIG. 2. The bottom plots show simulation results of the scaled average domain-wall velocity t" ' ' 1+<T ' vs. time for 
a = 0.5, 1.0, and 1.5 with open triangles, squares, and circles respectively. The upper plots show the scaled square-velocity, 
v 2 ^ t 2<T ^ 1+(7 \ vs. time with the filled triangles, open squares, and filled circles. Also shown, for comparison, is the t 1 ^ 3 line 
expected for the scaled square-velocity with a — 1.5. The systems have periodic boundary conditions, a fixed time-step, and 
start with 3200, 1600, and 800 domain walls for a = 0.5, 1.0, and 1.5 respectively. Note that the statistical errors shown are 
not independent. 



FIG. 3. 3>, the ratio of the measured to expected density (2.23), plotted vs. the number of domain walls (L^ /L(t)) on the 
left and vs. the ratio of measured to initial density (p(i)/p(0)) on the right. The data are for systems starting with 200, 800, 
3200, and 12800 domain walls for open triangles, squares, circles, and filled triangles, respectively. All of the systems have free 
boundary conditions and an initial random placement of domain walls. 



FIG. 4. The ratio of the measured to expected density up to a constant factor (2.23) plotted vs. the number of domain walls 
for a — 0.5. The top figure corresponds to periodic boundary conditions with no cutoff to the force law, while the bottom figure 



corresponds to free boundary conditions, 
starting with 1600 domain walls. 



The point types are the same as for Figure Bl with stars corresponding to systems 



FIG. 5. The ratio of the measured to expected energy density (2.24), up to a constant factor, plotted vs. the ratio of 
measured to initial density of domain walls. The top figure is for a = 0.5, while the bottom is for a — 1.0. The boundary 
conditions are periodic. The point-types are the same as the previous two figures. 



FIG. 6. The distribution of domain sizes, Fl(x) — Lrit(xL)/N, as a function of scaled domain size, x — l/L(t). Note that 
we have scaled with respect to the average domain size, L(t). The triangles, squares, and circles correspond to a — 0.5, 1.0, 
and 1.5, respectively. The filled points correspond to simulations with four times as many updates per factor of two in time. 
The l ine in the upper right is the exact result for a — ■> oo. The straight lines indicate power- law fits for small x (see equation 
(2.20)), as expected for a uniform rate of domain annihilation. 



FIG. 7. Real space correlations plotted against x = r/L(t). The top three figures are the scaled equal-time correlations 
f"(x) — L 2 C"(r,t) for a = 0.5, 1.0, and 1.5, while the bottom three figures are the corresponding scaled time-derivative 
correlations g(x) — (L/L) 2 T(r,t). The open circles correspond to correlations at times greater by a factor of 4 compares to 
triangles. The solid lines are described in the text: the upper curve represents a parameterless estimate of the correlations, while 
the offset straight line represents the expected asymptotic power-law for small x. The filled circles correspond to simulations 
with a finer timestep. The error bars have been suppressed for clarity. No data is shown for g(x) with x > 1/2. 



FIG. 8. For systems with short-range interactions, or with long-range interactions where a — > oo, the smallest domain will 
annihilate instantaneously with respect to larger domains. A single large domain is formed of size equal to the smallest domain 
and its two neighbours. The filled circles indicate domain walls. 



TABLE I. The scaling amplitudes of the grow th la w (se e 2.2S), where <7jj(<t) 
are for the dynamics specified in equations (2.2) and 



$» [2(1 + a 



|&. These amplitudes 



0.5 
0.75 
1.0 
1.5 



0.200 ± 0.005 
0.180 ±0.005 
0.175 ±0.005 
0.170 ± 0.005 



0.42 ±0.01 
0.37 ±0.01 
0.35 ±0.01 
0.32 ±0.01 



13 



oo 



e -7E/ 2 ~ 0.28 



g-7B/ 2 ~ 0.28 



TABLE II. The coefficient of the velocity of fast-moving domain walls, v(l) = j3l " for I <C L, obtained by fitting figure 
by the distribution in equation (|2.19|). For a — > oo, the domains are uncorrelated and /? = 1 (see section 



a 


13(a) 


0.5 


0.75 ± 0.05 


1.0 


0.90 ± 0.05 


1.5 


0.95 ± 0.05 


oo 


1 
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